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Abstract A new method for automated detection of polar coronal holes is pre- 
sented. This method, called perimeter tracing, uses a series of 171, 195, and 304 A 
full disk images from the Extreme ultraviolet Imaging Telescope (EIT) on SOHO 
over solar cycle 23 to measure the perimeter of polar coronal holes as they appear 
on the limbs. Perimeter tracing minimizes line-of-sight obscurations caused by 
the emitting plasma of the various wavelengths by taking measurements at the 
solar limb. Perimeter tracing also allows for the polar rotation period to emerge 
organically from the data as 33 days. We have called this the Harvey rotation 
rate and count Haxvey rotations starting 4 January 1900. From the measured 
perimeter, we are then able to fit a curve to the data and derive an area within 
the line of best fit. We observe the area of the northern polar hole area in 1996, 
at the beginning of solar cycle 23, to be about 4.2% of the total solar surface 
area and about 3.6% in 2007. The area of the southern polar hole is observed 
to be about 4.0% in 1996 and about 3.4% in 2007. Thus, both the north and 
south polar hole areas are no more than 15% smaller now than they were at the 
beginning of cycle 23. This compares to the polar magnetic field measured to be 
about 40% less now than it was a cycle ago. 

Keywords: Coronal Holes; Solar Cycle, Observations 
1. Introduction 

The Sun's polar magnetic field has been used to forec;ast upcoming solar activity 
for four solar cycles (Schatten et ai, 1978; Schatten, 2005; Svalgaard, Cliver, 
and Kamide, 2005). The solar magnetic field that extends far beyond the corona 
into the interplanetary medium is related to the area of the corona threaded 
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by the field. These regions are called solar coronal holes and are defined by 
both their low emission in EUV and soft x-ray wavelengths and density in the 
corona as well as their dominantly unipolar magnetic fields (Bohlin, 1977). First 
observed in 1941 above the solar limb in visible light (Waldmeier, 1981), coronal 
holes are most easily seen in space-based EUV or x-ray images. For ground- 
based measurements, coronal holes are visible in the equivalent width of the 
He I 10830 A line (Andretta and Jones, 1997). 

Polar coronal holes are coronal holes that occur at the north and south pole of 
the sun and appear as dark caps on EUV images. Polar holes deserve particular 
attention because, as the solar cycle progresses toward sunspot number {Rz) 
maximum, the polar coronal holes become poorly defined and ultimately disap- 
pear. During the decline from solar maximum, the polar coronal holes (PCHs) 
reemerge, becoming better defined and almost symmetric at solar minimum. 
Due to their relatively long life, PCHs are the source of the majority the solar 
interplanetary quiet magnetic field. 

On average, the polar magnetic field behaves in much the same way as the 
PCHs. The polar field strength at both poles reaches a maximum at solar min- 
imum. At solar maximum, the polar field decreases through zero on its way to 
changing its polarity. Over the past three solar minima the magnitude of the 
polar field has steadily decreased, indicating a decreasing trend in solar activity. 
The current solar minimum has a polar field magnitude of 0.6 gauss, compared 
to 1.0 gauss in the minimum preceding Cycle 23 as shown in Figure 1 (Hoeksema, 
2008). 

Coronal holes, both polar and equatorial, have historically been identified 
visually and hand-traced by experienced observers. Examples of these images 
can be seen in Harvey and Recely (2002) and at the National Solar Observatory 
(NSO) Synoptic Optical Long-term Investigations of the Sun (SOLIS) project. 
More recently, there have been a few attempts to automate the process of coronal 
hole identification and detection (de Toma and Arge, 2005; Henney and Harvey, 
2005; Malanushenko and Jones, 2005; de Wit, 2006). 

Both manual and automatic methods in detecting PCHs have difficulties. 
Manual dctccition methods are extremely time consuming and difficult to quan- 
tify, as they rely on the observers' qualitative assessments of the images. Many 
automated methods use synoptic maps for hole detection (Harvey and Recely, 
2002). However, synoptic maps do not lend themselves to polar analysis. Syn- 
optic maps are compiled from sections of images of the solar disk, taken over a 
Carrington rotation (CR) of 27.27 days, and then mapped onto a rectangular 
grid (Benevolenskaya, Kosovichev, and Scherrer, 2001). Projecting a sphere into 
a 2D map is always problematic due to the resulting distortions in the apparent 
size of surface features; (e.g., the size of Greenland on a Mercator projection of 
the Earth). On the sun this problem is compounded by solar differential rotation. 
The polar regions of the sun rotate once every approximately 33 days. Using syn- 
optic maps to investigate polar features forces a CR period onto those features. 
Due to these problems, a better method to measure PCHs was developed, as 
discussed below. 
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2. A New Way of Looking at PCHs 

The appearance of coronal holes differs in each wavelength, which makes it 
unreliable to use one wavelength alone to measure them. Solar and Heliospheric 
Observatory's (SOHO) Extreme ultraviolet Imaging Telescope (EIT) observa- 
tions over the last solar cycle gives us a unique opportunity to have extreme 
ultraviolet images of the sun in four different wavelengths at nearly the same 
moment in time (Moses et al., 1997). EUV Fe lines measured by EIT have 
good contrast, but the emitting material has a relatively large scale height; for 
example, 195 A Fe xii has an emission that peaks at a temperature of 1.6 x 10^ 
K and has an associated scale height between 60 and 190 Mm (Moses et al., 
1997). When looking at the polar region on the solar sphere, if the emitting 
plasma has a small scale height, the holes are more easily discerned. However if 
the emitting plasma has a relatively large scale height, the edge of the polar hole 
is more easily seen at the limb of the the sun. In effect, the large scale-height 
is obscuring the edge of the polar hole at the central meridian (Figure 2). This 
would naturally lead one to choose a wavelength with the smallest scale height 
to get the most accurate measurement of the perimeter of the PCH. The He II 
304 A line has a small scale height compared to the other lines that EIT measures 
(He II 304 A has an emission that peaks at a temperature of 8.0 x 10^ K and 
has scale height between 3 and 10 Mm). However, the corrugations in the 304 A 
chromospheric network have a similar contrast to coronal holes, which makes 
discerning the boundary of the holes from the quiet sun network difficult. The 
high-temperature Fe xv 284 A images from EIT are particularly problematic 
to analyze, because at solar minimum, hot plasma is absent, resulting in poorly 
defined 284 A images (Kahler and Hudson, 2002). Thus, despite their limitations, 
the Fe lines at 171 and 195 A and the He ii line at 304 A images are the best 
EUV images to use (Figure 3). The He ii 10830 A images from Kitt Peak could 
be used to supplement the EIT images; however they are left out of this survey 
to simplify the discussion. 

Because of the differential rotation rates of the sun, it would be inaccurate 
to use the Carrington rotation rate as the polar rotation rate. The mean solar 
rotation rate, co (in degrees per day) , is a function of heliographic latitude: 



where A, B and C are empirically determined coefficients and <j) is the latitude. 
The values of the coefficients used are listed in Table 1 (Howard, 2000) . Coronal 
features, surface plasma and magnetic field features were used in determining the 
rotation rates of PCHs, because properties of coronal holes can fit into all three. 
The average rotation rate of these features is 10.88° per day or just slightly over 
33 days for one rotation at a latitude of 70°. Thus we approximated the PCH 
rotation rate to be 33 days, from here forward called the Harvey rotation rate. 
For this work, we have defined the Harvey Rotation number (HR) to be 



u) = A + Bsm^ (j) + Csin* (j), 



(1) 



HR = 



JD - 2415024.5 
33 



(2) 
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where JD is the JuHan day and Harvey rotation zero began 4 January 1900 (the 
beginning of Carrington Rotation 619). Similarly, hcliographic coordinates can 
be translated into an HR adjusted coordinate system by letting 0° in the HR 
coordinate system be the central meridian of the sun at the beginning of the 
Harvey rotation. 

3. Detection Method 

To determine the size of the polar coronal holes on sun using perimeter tracing, 
while avoiding some of the obscniration of the hole boundary when working in 
the polar region, the method of Henney and Harvey (2005) is adapted to a limb 
geometry. One image per day, as close to 0:00 UT as possible, is selected for each 
of the three EUV wavelengths. 

3.1. Calibration 

The images are first calibrated with the eit_prep software using the calibration 
tables available in the Solarsoft suite (Figure 3). They are also truncated so 
that the raw data numbers are whole numbers that fall between and 1000. 
Because coronal holes have low emission in the EUV, this truncation does not 
affect the detection of the hole but does reduce the storage requirements of the 
image array. 

3.2. Image Analysis 

The image is then spatially smoothed. The morphological image analysis (MIA) 
functions Dilate (V) and Erode {£) are applied in series to remove small noisy 
struc;tures. In general, MIA functions can be thought of as a 2D convolution of 
the shape operator and the original image. Erode and Dilate are defined as: 

V{I, S) {i, j) = max -k,j-l) + S{k,l)] (3) 

<k,l>s 

and 

£{I,S){t,j)= min [I{i + k,j + l)-S{k,l)] (4) 

<k,l>s 

where is the original image array and S{k,l) is a shape-operator array. 

The shape-operator array is defined by its displacement, <k, l>s, from its origin, 
(0,0) (Michielsen and Raedt, 2001). 

Erode and Dilate can be combined to create two other MIA functions Open 
(O) and Close (C), which are defined as: 

0{I,S){i,j) = V[£{I,S),S] (5) 

and 

C{I,S){i,j) = £[V{I,S),S]. (6) 
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The functions Open and Close are also idempotent, 0{0{I,S),S) = 0{I,S) 
and C{C{I, S), S) — C{I, S), a property that is advantageous in scientific image 
analysis. Generally, the Close function fills in small gaps between large regions 
of the same intensity. This has the effect of smoothing over regions that are too 
small to be part of a coronal hole while leaving larger regions largely intact. 
Features that do not fit the shape-operator (5) passed through the image are 
removed. The Open function, on the other hand, associates small features in 
close proximity to each other effectively grouping small regions of the same 
intensity into one large region. On solar images, using the Open function creates 
smooth edges between features. For a further discussion of MIA functions and 
their properties, see Michielsen and Raedt (2001). 

In this c;as(;, we used a filter (J-') having a constant shape-operator [S) of a 
circular kernel with a radius of 3 pixels. The Filter is defined by: 

J^{I,S){i,j)=C[0{I,S),S], (7) 

and has the effect of smoothing out the subtle variations in the image while 
leaving larger image features intact (Figure 4A). These image manipulations are 
needed to eliminate bright points that can occur in coronal holes and distract 
the hole boundary detection. It is also needed to blur out the grain structure in 
304 A images and reveal the hole boundary. 

3.3. Thresholding 

The sun's intensity changes throughout the solar cycle. When detecting PCHs, 

this can be a problem, because the same intensity value may be a maximum 
in one image and a minimum in a different image of the same wavelength. To 
overcome this problem, an integrated intensity histogram of the remaining image 
is created. The integrated histogram is then divided by the total number of pixels 
in the image. From this histogram, a threshold level is picked to reveal the 
candidate holes in the disk (Figure 4B) . The threshold level is different for each 
wavelength. For example, in a 171 A image, the threshold is the intensity-bin 
whose value is 0.70, or where 30% of the pixels are darker than that intensity- 
bin. In a 195 A image, the bin value is 0.73; and in a 304 A image it is 0.50. 
These bin values were derived empirically to best expose the polar holes in each 
of the wavelengths. A binary image is then produced using the threshold level 
determined by the integrated intensity histograms where values less than the 
threshold is set to zero and all others are set to 1. 

Notice that, in the binary threshold image produced, exposing the polar holes 
on the limb of the sun overexposes holes near the disk center. The thresholding 
values distinguish this process from an equatorial detection method, because 
we set the threshold to specifically fit our limb geometry. Since we are only 
concerned with the limb, overexposing the equatorial holes is not a concern. By 
changing the threshold values we could easily adapt this algorithm to detect 
coronal holes in other regions. 
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3.4. Extraction 



Next, wc isolate the disk itself by trimming the coronal streamers outside the 
limb of the sun. The location of the solar limb is determined using limb bright- 
ening as defined by the EIT package in SolarSoft. The center of the disk is also 
trimmed, leaving an annulus spanning the outer 6% of the solar radius. Then 
the polar regions are isolated by masking the disk up to 60° north and south 
latitude (Figure 4C). The location of the PCH boundary at the limb can now 
be extracted. 

Once the candidate hole edges are identified, their heliographic longitudes are 
calculated and translated into Harvey coordinates (Figure 4D). Using a Harvey 
coordinate system to determine the edge of the holes has the added benefit of 
allowing the perimeter of the hole to emerge organically from the data; i.e. the 
perimeter tracks back onto itself after a 33 day rotation. The points are then 
plotted on a polar grid using the HR rate (Figure 5) . 

To evaluate the quality of the measurement, we find the fraction of pixels in 
the remaining binary limb within the candidate hole edge that appear as a hole 
(Pixelszeio) as compared to all of the pixels (Pixelstotai) that fall within that 
candidate hole boundary, i.e., within the boundaries identified in Figure 4D, the 
number of white pixels is compared to the total pixels in Figure 4C (eq. 8). 



Quality = -r^-^ (8) 
Prxelstotai 

If the hole is poorly defined, then the quality is low. If the quality of the 
measurement of the hole boundary is less than 15% then the measurements 
are discarded. This is to help guard against false positive measurements due to 
equatorial holes encroaching into the polar region, especially when the sun is 
active. 



3.5. Method Verification 



To test the perimeter tracing algorithm we generated a pseudo-coronal hole on 
a computer-generated sun. The circular hole was 20° offset from the rotational 
axis, 18° in diameter beveled up to the surface (the hole is about 23° in diameter 
at the surface), and the rotational axis was inclined -7° from our visual vertical 
axis (simulating the most extreme heliographic tilt seen during the year. Fig- 
ure 6A). We rotated the model through 360° in 15° steps letting the detection 
program find the edges of the simulated hole every step. Then, after assigning 
spherical coordinates to each of the detected points, we mapped them to a polar 
grid (Figure 6B). The result returned a hole centered at a colatitude of 23° 
with a diameter of 20°. The diameter is not exactly 18° because the detection 
method is sensitive to the bevel of the simulated hole. Also, the detected limb 
of the simulated hole is not exactly circular because of the projection effects as 
the hole sweeps around. 
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4. Data Analysis 

When looking at an HR of polar measurements, it would be imprudent, although 
tempting, to simply connect the data points and say that this represents a good 
measurement of the PCH boundary; however, it would tend to emphasize the 
effects of the noise in the measured data and would produce a less accurate 
area bounded by the limb measurements of the hole perimeter. In Section 3 we 
accounted for this by using morphological transforms to smooth the candidate 
hole boundaries. When assessing the area within the coronal hole boundaries, 
fitting the collected data with a line of best fit keeps the integrity of the physical 
nature of coronal holes. 

4.1. Functional Fitting 

The first step in the fitting process is to combine the two separate limb measure- 
ments into one data set. This is where translating into Harvey coordinates in 
Section 3.4 pays off. Combining the sets of points is as simple as creating a large 
array from the two smaller arrays, since Harvey coordinates are the appropriate 
reference frame for polar features. 

A first order filter is then applied. If greater than half of the measured points 
are missing in a given Harvey rotation, the entire rotation is thrown out. This 
filter is needed to guard against under-sampling the limb measurements of the 



The array is then translated so that the center of mass of the points is the 
origin of the coordinate grid. This translation is necessary to force the fitting 
routine to fit each point more efficiently. For example, if the polar hole intersected 
the origin, the fitting routine would tend to extend into negative space to fit the 
feature, creating an unrealistic fit of the hole. Translating the PCH measurements 
provides some protection against this circumstance. 

The array of PCH measurements can then be fitted with a Fourier cosine 
series in the longitude: 



-^max is allowed to vary from two to seven to generate several candidate fits. 

The points are fitted using a non-linear least squares fitting routine with an 
arbitrary number of input parameters. This method uses an iterative gradient- 
expansion algorithm combined with a minimization of (Press et al., 1996). A 
full discussion of the fitting procedure can be found in Bevington and Robinson 
(1992). A similar method was used by (Pesnell et al., 2000) to analyze ozone 
variations in the Earth's atmosphere. 

To increase the viability of the fit produced by the procedure described above, 
weights are assigned to each measurement. Because PCHs do not have disconti- 
nuities, we can say adjacent measurements are more accurate if they do not have 
large variations. In order to give greater weight to points that are more closely 
associated while putting less emphasis on outlying points, a Gaussian weighting 



PCH. 




(9) 
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system is used. The total array is divided up into octants: between 0° and 45°, 
45° and 90°, etc. The variance, cr^, is then calculated within each octant, for 
the angular distance to the pole {9) where the standard deviation in the radius 
is (Bevington and Robinson, 1992): 



Since points that have a small radial variance mean that the edge of the PCH is 

better defined, wc can weight the points as the reciprocal of the the variance in 
each octant. Note that all the points in an octant are given the same weighting. 

Since we generate six different lines of best fit using Equation (9), we must 
make some judgments about which fits arc better than others. If the fitting line 
(/[0]) extends below 45° latitude, the line is no longer fitting in the polar region 
and we can disregard that fit. Next, we must be careful of over-fitting the given 
polar measurements. A symptom of over-fitting is a large length of the curve 
(i.e., it has a lot of wiggles) that defines the perimeter of the PCH. The arc 
length of the fitted curve is found analytically by integrating the generated line 
of best fit (/[</>]). If a fit has a total length that is greater than four times the 
smallest curve's length, we disregard the larger fit. 

Of the remaining lines of best fit for the given PCH data, we can now find 
the area. The area of the polygon formed by the line of best fit can be found by 
integrating /(0) from Equation (9) (Beyer, 1991). Another affect of over-fitting 
a series of points is that a larger than necessary area is enclosed. Thus we can 
eliminate all fits that are greater than twice the smallest enclosed area. 

From the fits that remain, taking a simple average of the areas they enclose 
gives us a good sense of the fitted area of the PCH (the black line in Figure 5). 

4.2. Statistical Analysis 

It is important to define a confidence measurement when working with fits to 
a set of data. To accomplish this, we estimate the probable uncertainties in 
the fitted parameters. Because each of the data points is independent, each 
contributes its own bit of uncertainty to the fitting parameters. Propagating 
these errors, it can be shown that the variance a'j in the value of the function 



Directly evaluating the partials of each of the fitting parameters from the solution 
will produce the variances in each of the parameters (Press et al, 1996). 

With the variances in each parameter, we can then produce a one-sigma 
confidence range. By subtracting one sigma from each of the parameters, a 
minimum area fit is produced. Likewise, by adding one sigma, a maximum area fit 
can also be determined. However, since the fitting function is a transcendental 
function, the error range does not always produce realistic results. When the 




(10) 



will be 




(11) 
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fitting function has a tenuous relationship with a few outlying points, adding or 
subtracting a onc-sigma error can cause significant ringing in the fitting function 
leading to an unrealistic fit. This problem is most evident during solar maximum 
where the measured points are not well correlated. 

To determine the range of values seen in Figure 7, we found the maximum 
and minimum areas associated with each of the fits determined to be good in 
Section 4.1. Using the same method as used with finding the area bounded 
by a line of best fit, we averaged the maximum values and minimum values 
respectively. This produced the one-sigma range of values. Where the errors 
were not realistic, the area bounded by the best fit was plotted. 

5. Results and Discussion 

Studying the evolution of the polar coronal holes over solar cycle 23 utilizes 
the entire lifetime of SOHO and chronicles tiic (kH;line of the PCH into solar 
maximum and the rise to the current solar minimum (Figure 7). The northern 
PCH disappears in the final months of 1999 and returns in the last quarter of 
2000. On the other hand, the southern PCH does not disappear until mid 2000 
and does not return in earnest until the first quarter of 2002. This asymmetry 
is consistent with past surveys of PCHs. 

The northern PCH is significantly better defined than the southern hole near 
solar maximum, as seen by the large error in the southern measurements in 2002 
and 2003. The large error measurements in the south are confirmed by a qual- 
itative assessment of the disk images during 2002 and 2003. These disk images 
appear to have large equatorial holes migrating toward the southern polar region 
making it difficult to quantify the boundary of the PCH. It is also interesting 
to note the relative disagreement of PCH area of the three wavelengths in 1996 
and into 1997 in contrast to their relative agreement (within each others error 
estimates) for the rest of the series. 

Because of the relatively small time series of this study, it is difficult to make 
broad assertions about the nature of the polar holes. Nonetheless, comparing 
this series to previously compiled data reveals some notewortiiy results. 

Comparing this work to that of Harvey and Recely (2002) reveals some strik- 
ing similarities. Harvey and Recely produced a time-series plot of the PCH area 
over the second half of solar c;ycle 22 and first half of 23. Harvey and Recely 
have the northern polar hole shrinking to less than 1% of the surface area of 
the sun midway through 1999 and Harvey's "pre-polar hole" reemerging at the 
beginning of 2001. They also show the southern PCH disappearing completely 
in mid 2000 but end their study before the southern hole reappears. Figure 7 
shows a similar time evolution of the northern and southern holes. 

After the northern PCH reemerges, Harvey and Recely (2002) observed an 
area of about 4.8% of the surface area of the sun. By comparison, this work 
observes a northern area of about 3.5% or an He I 10830 A area that is about 
37% greater than the EUV-derived area. In 1996, Harvey and Recely report the 
northern PCH area to be about 8.25% the solar surface area. We observe an 
area of about 4.2%, meaning that the He I 10830 A area is almost twice the 
EUV-derived area at the beginning of solar cycle 23. 
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Another important comparison is with the Wilcox Solar Observatory (WSO) 
measurements of the polar magnetic field (Figure 7). The WSO reports a mag- 
netic flux magnitude of about 1.0 gauss at the minimum preceding cycle 23, 
while it reports about 0.6 gauss at the current minimum. The polar field is 40% 
less now than it was a cycle ago. Because coronal holes arc characteristically 
threaded by open magnetic field lines and the solar polar magnetic fields are 
confined to relatively small caps (Sheeley, Wang, and Harvey, 1989), this seems 
imply that the size of the polar hole should be significantly smaller as well. 
However, we observe the size of the hole to be at most 15% smaller this solar 
minimum as compared to last. This discrepancy would suggest that the polar 
coronal hole is more a product of dynamics rather than magnetic fields. 

When comparing the WSO time series of the magnetic field to the area of the 
polar holes at both poles, it appears that there is a close correlation between 
them. A correlation plot of the north or south hole and the WSO field reveals 
no more than a weak correlation between the plots seen in Figure 7. If we 
separate the average WSO field into its northern and southern components, 
the correlation seen is about the same as with the average magnetic field. 

6. Comments and Conclusions 

We have shown that perimeter tracing on disk images liberates the analysis of 
PCHs from a reliance on synoptic maps. Perimeter tracing in the polar region, 
given the differential coronal rotation rate, demands the development of a new 
coordinate system rather than a mid-latitude Carrington rotational period. A 
Harvey coordinate system better represents the polar rotation rate and allows 
us to get a more accurate representation of the PCH. Using EIT images over 
solar cycle 23, we observe the PCH area of the north at the end of 2007 to be 
about 3.6% of the total solar surface area versus about 4.2% in 1996. The area of 
the southern PCH at the end of 2007 is observed to be about 3.4% as compared 
to about 4.0% in 1996. The estimate of the PCH size between 1996 and 1999 is 
complicated, because the three EIT wavelengths used have distinctly different 
areas. After the beginning of 1999 the areas of the different wavelengths agree 
more closely. Yet even with this difference, we observe both the north and south 
PCH area to be no more than 15% smaller now than they were at the beginning 
of cycle 23. 

It is possible to incorporate magnetograms into this PCH detection method 
by summing the net magnetic field within the derived hole boundary. Using 
magnetograms would allow us to measure the magnetic fltix density evolution 
over the solar cycle and better isolate the polar hole from equatorial holes. 
Perimeter tracing is also extendable into the past with archives of He I 10830 A 
images, which will allow us to get a better idea of the long-term evolution of 
PCHs. 
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Figure 1. The average magnetic field measured at WSO beginning in July of 1976 through 
the end of 2007. Measurements are made every 10 to 14 days of the field at both poles. This 
plot shows the average of these measurements from both poles, (Sn — Bs)/2, smoothed with a 
7 point boxcar function to remove the effects of the tilt of the solar rotation axis with respect 
to the ecliptic. The times when the averaged magnetic field is zero correspond to the polarity 
fiip associated with solar maximum. 
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Figure 2. A diagram showing the effects of the emitting plasma's scale height on the visibility 
of the edge of the polar coronal hole. The view shows a slice of the northern hemisphere of 
the sun at the central meridian; the rotation axis is also marked. A) shows an example of a 
relatively small scale height contrasted with a large scale height of B). Notice how the emitting 
plasma in B) obscures the edge of the polar hole. 
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Figure 3. Examples of calibrated SOHO/EIT images at different wavelengths where polar 
coronal holes (PCHs) can be seen. A) 171 A, B) 195 A, and C) 304 A were all taken near 
June 1, 2007. The PCHs in these wavelengths appear as dark polar caps in both the north and 
south. The sizes and shapes of the PCHs vary between each wavelength due to differing scale 
heights in the emitting plasma. 
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Figure 4. An illustration of the automated PCH detection process using the 195 A image 
shown in Figure 3: A) After using morphological transform functions Close and Open to blur 
the original image; B) A binary image of the image shown in A) retaining the top 73% of the 
bin values in an integrated intensity histogram; C) An annulus of the outter 6% of solar disk 
after removing off-limb and central disk data; D) The edges of the north and south PCHs are 
shown with a circle. The heliographic coordinates are then calculated and used to mark the 
perimeter of the PCH after translation to a Harvey coordinate system. 
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Figure 5. The north and south pole limb measurements in Harvey coordinates beginning 12 
January 2006. Each blue and red mark designates the coordinates of the edge of the PCH 
on the ascending and descending limb respectively. A fit to the data using Equation (9) with 
A^max = 7 is drawn as the solid line. The center of the PCHs is coincident with the rotation 
axis. 





Figure 6. Example of converting a simulated PCH to a polar projection using limb detection 
and perimeter tracing. A) One of the images of the simulated sun with a northern PCH. The 
simulated hole has an offset of 20° from the pole and a diameter of 18° at its lowest point. The 
whole simulated sun has a heliographic tilt of -7°. B) The perimeter found by the perimeter 
tracing technique described in Section 3. The PCH at the ascending (western - blue) and 
descending (eastern - red) limb of the sphere are shown. The detected hole has a diameter of 
20°, and its center (marked with a star) has an offset of 23°. 
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Figure 7. A time series of PCH areas as a percent of the total solar surface area. Each area is 
calculated over a 33-day Harvey rotation using the perimeter tracing technique. The blue line 
corresponds to 171 A, the green to 195 A, and the red to 304 A. All lines are smoothed with 
a 3 point boxcar function. The yellow hash-marked sections on the graphs mark times when 
more than half of the Level EIT data was missing for a given Harvey rotation. The absolute 
value of the averaged WSO polar magnetic field from Figure 1 is included for reference. 



SOLA: PCH_Draft_v3.tex; 23 January 2009; 14:44; p. 17 



Kirk ct al. 



Table 1. The coefficients used in Equation (1) and the 
resulting rotation rate (cd) for each type of feature to 
determine the polar rotation rate of PCHs (Howard, 2000) . 

Feature ABC u> 

deg day~^ 

Coronal features 13.46 -2.99 — 10.82 
Surface plasma 14.11 -1.70 -2.35 10.77 
Magnetic field 14.37 -2.30 -1.62 11.07 
Average 10.88 
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